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^ ■ Abstract 

CO ' We study the low- momentum behaviour of Yang-Mills propagators obtained from Landau- 

gauge Dyson- Schwinger equations (DSE) in the PT-BFM scheme. We compare the ghost 
propagator numerical results with the analytical ones obtained by analyzing the low-momentum 
f~| . behaviour of the ghost propagator DSE in Landau gauge, assuming for the truncation a con- 

On' stant ghost-gluon vertex and a simple model for a massive gluon propagator. The asymptotic 

("^ . expression obtained for the regular or decoupling ghost dressing function up to the order 

g_) ' C'(5^) is proven to fit pretty well the numerical PT-BFM results. Furthermore, when the 

r^ . size of the coupling renormalized at some scale approaches some critical value, the numerical 

PT-BFM propagators tend to behave as the scaling ones. We also show that the scaling solu- 
psj , tion, implying a diverging ghost dressing function, cannot be a DSE solution in the PT-BFM 

^ ' scheme but an unattainable limiting case. 
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^ ! 1 Introduction 

O 



The low-momentum behaviour of the Yang-Mills propagators derived either from the tower 
of Dyson-Schwinger equations (DSE) or from Lattice simulations in Landau gauge has been a 
very interesting and hot topic for the last few years. It seems by now well established that, if we 
5^ I assume in the vanishing momentum limit a ghost dressing function behaving as F^q^) ~ (f7^)"^ 

and a gluon propagator as A(g^) ~ {q'^)'^'^~^ (or, by following a notation commonly used, a 
gluon dressing function as G{q^) = (f'/S.{(f') ~ {(f')'^^), two classes of solutions may emerge (see, 
for instance, the discussion of refs. [Il[2]) from the DSE: (i) those, dubbed "decoupling", where 
ap = and the suppression of the ghost contribution to the gluon propagator DSE results in 
a massive gluon propagator (see [3ll3j and references therein); and (ii) those, dubbed "scaling", 
where ap 7^ and the low- momentum behaviour of both gluon and ghost propagators are related 
by the coupled system of DSE through the condition 2aF + o-g = implying that F'^{q'^)G{q^) 
goes to a non- vanishing constant when g^ — )> (see [SJE] and references therein) . As a matter of 
fact, F^(g^)G(g^), which gives the perturbative running for the coupling constant renormalized 
in Taylor-scheme [7], is a good quantity in discriminating the kind of solutions we deal with. It 
is worth to remember that, despite the widely accepted nomeclature for the classes of solutions, 
neither a scale invariance nor a decoupling of the IR dynamics for the theory can be inferred 
from the low-momentum behaviour of such a Taylor coupling. 



How both types of IR solutions for Landau gauge DSE emerge and how the transition between 
them occurs, being governed by the size of the couphng taken as an integration boundary condi- 
tion at the renormalization momentum, was initially discussed in ref. [1] through the analysis of 
a ghost propagator DSE combined with a gluon propagator taken from lattice computations. It 
should be remembered that one needs to know the QCD mass scale to predict the QCD coupling 
at any momentum. This mass scale should be of course supplied to get a particular solution 
from DSE and can be univocally related to the boundary condition needed, after applying a 
trunctation scheme, to solve the equation. The existence of a critical value for the coupling at 
any renormalization momentum was suggested by that partial analysis. No solution was proved 
to exist for any coupling bigger than the critical one and the unique scaling solutiorQj appeared to 
emerge when the coupling took that critical value. Later, the authors of ref. |6] confirmed, by the 
analysis of the tower of DSE truncated within two different schemes and also in the framework 
of the functional renormalization group, that the boundary condition for the DSE integration 
determined whether a decoupling or the scaling solution occurs. A similar analysis have been 
recently done in the Coulomb gauge [8] leading to the same pattern of ref. [1], although the 
authors interpreted the boundary condition in terms of the gauge- fixing ambiguity (see also [9] ) . 
Furthermore, an analytic study based on the pinch technique in ref. [10] shows that, within some 
approximations, there is a lower limit for the gluon mass, which the IR singularities of QCD 
still persist below, that can be also interpreted as an upper limit to the coupling. Very recently 
also, a next-to-leading low-momentum asymptotic formula for the decoupling ghost dressing 
function solutions was obtained by studying the ghost propagator DSE with the assumption, 
for the truncation, of a constant ghost-gluon vertex and of a simple model for a massive gluon 
propagator [TT]. In this asymptotic formula, the ghost-propagator low- momentum behaviour 
appeared to be regulated by the zero-momentum effective charge in Taylor scheme [12] and by 
the Landau-gauge gluon mass scale. 

In the present note, the work of ref. [1] will be extended by the analysis of the results [13] 
obtained by solving the coupled system of Landau gauge ghost and gluon propagators DSE 
within the framework of the pinching technique in the background field method [13] (PT-BFM) . 
Our main goal is investigating whether the same pattern for regular (decoupling) and critical 
(scaling) solutions in ref. [1] is also found for the PT-BFM solutions. In the PT-BFM scheme, the 
zero-momentum ghost dressing function will be seen to diverge when the coupling approaches 
some critical value, as it should be expected for the scaling solution {ap 7^ 0). This seems to 
support the suggestion of a transition from one to another solutions controled by the size of 
the coupling approaching a critical value [T]. The authors of ref. [S] obtained similar results 
but they applied the zero-momentum ghost propagator as the boundary condition for the DSEs 
integration and missed its connection with the value of the coupling at the renormalization 
momentum (i.e. the particular value of Aqcd one applies to build the solutions) or the critical 
coupling the scaling behaviour requires to emergqj. This connection is an important ingredient 
because it provides us with a manner, through a comparison with the physical strong coupling, 
to discuss whether the scaling critical DSE solution could be allowed by the data. We will 
also compare the low-momentum analytic results of ref. [11] with the PT-BFM results. In 
particular, the next-to-leading asymptotic formula for the ghost dressing function with ap = 



^The authors of [TS] proved there, once the scahng behaviour is assumed, the uniqueness for Yang-Mills infrared 
solutions 

^These authors furthermore invoked the renormalization group invariance to claim that such a critical coupling 
does not appear. Nevertheless, the renormalization group invariance only requires for the critical values of the 
coupling at any two fixed renormalized momenta to be connected by the appropriate renormalization group 
running (the same, of course, for any decoupling solution). Indeed, the renormalization flow for the coupling can 
be defined by a{q^) = a{iJ?)F^ {q^)G{q^) to satisfy this and it is found to agree, in the perturbative domain, with 
the perturbative running given by the /3-function in Taylor scheme [7]. 



will be shown to nicely describe the low-momentum PT-BFM results for different values of the 
coupling taken at the renormalization point, as a boundary condition for the DSE integration. 
It should be emphasized that the ansatz for the massive gluon propagator applied to derive 
the analytical results in ref. [11] has nothing to do either with the numerical analysis based on 
lattice results in [1] or with the one based on PT-BFM results in this paper, both leading to 
obtain a critical coupling. These analytical low momentum results are applied in this paper only 
for comparative purposes. However, as a result of the comparison, this boundary condition will 
be put in direct relation to the zero-momentum effective charge in Taylor scheme |12j and to 
the Landau gauge gluon mass. Finally, we will also argue that a scaling solution cannot exist 
as a solution of the coupled system of PT-BFM DSE, although the PT-BFM solutions tend to 
it when the coupling aproaches the critical value. More generally, a diverging ghost dressing 
function cannot be obtained when the gluon propagator is massive [ac = 1)- Of course, this 
last claim is not a new result: it is well-known that the scaling solution only can emerge if 
the ghost dominance in the gluon propagator DSE, after assuming ap < 0, implies the gluon 
propagator to vanish at zero momentum [5l[6]. However, it is not worthless to emphasize that, 
for the gluon propagator, being massive implies not to observe the relation lap + ac = and is 
therefore a sufficient condition for the scaling solution not to appear. It should be also pointed 
out that LQCD results (see |16H21] and references therein), pinching technique results (see, for 
instance, |10 ^ I14 ^ [22 | I23]). refined Gribov-Zwanziger |f| formalism (see [28]) or other approaches 
like the infrared mapping of A0 and Yang-Mills theories in ref. [29j or the massive extension of 
the Fadeev-Popov action in ref. [3D] appear to support a massive gluon propagator. 

We organized this note as follows: first, we briefly review the low-momentum behaviour for 
the propagator solutions of the ghost propagator DSE in section [21 we then compare the PT- 
BFM results with the low-momentum analytical expression and discuss their dependence with 
the size of the coupling at the renormalization point, taken as a boundary condition for DSE 
integration, in section [Sj and we finally conclude in section [H 

2 The two kinds of solutions of the ghost propagator Dyson- 
Schwinger equation 

As was explained in detail in refs. [2l[Tl], the low-momentum behavior for the solutions 
of the Dyson-Schwinger equation for the ghost propagator (GPDSE), which can be written 
diagrammatically as 

q-k 
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can be obtained in Landau gauge by procceeding as follows: 

we consider eq. ([T]) for two different (although parallel) external ghost momenta, k and p, 
such that p^ — k"^ = J^fc^ {5 being an extra parameter that, for the sake of simplicity, will be 
taken to be small enough as to expand on it around 0) and subtract them, as a regularization 
prescription for not to have to deal with any UV cut-off. Then, the ghost dressing function, 
F{k'^), and the gluon propagator form factor, A(A:^), are renormalized by applying the MOM 
prescription. 



•^In addition, K-I. Kondo triggered very recently an interesting discussion about the Gribov horizon condition 
and its implications on the Landau-gauge Yang- Mills infrared solutions [24H27| . 
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where /i^ is the subtraction point; and, as explained in [21[TI], we choose for the ghost-gluon 
vertex, 
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to apply the MOM prescription in Taylor kinematics {i.e. with a vanishing incoming ghost 
momentum) and assume the non-renormalizable bare ghost-gluon form factor, Hi{q,k) = Hi, 
to be constant (at least in the low-momentum regime for the incoming ghost). Thus, after being 
cast into a renormalized form and the above-mentioned subtraction, the GPDSE reads 
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and where the ghost dressing function and gluon propagator should be understood as renor- 
malized at the subtraction momentum, /i^, and 5r(/U^) is the gauge coupling renormalized at a 
subtraction point with Taylor kinematics (Taylor MOM scheme). Then, we cut integral domain 
of I^k"^) in eq. ([5]) into two pieces by introducing some new momentum scale, ^q, below which 
both ghost and gluon are assumed to be well described by the following ansatze 
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Thus, we shall look for the ghost dressing function, Fir, its leading behaviour being parameter- 
ized through a general power law behaviour where ap > —2 to keep the integral I{k'^) infrared 
convergent, and for a massive u gluon propagator, that implies of course a power law with 
ac = 1, as the current lattice data seems to point to. It is well-known that the low-momentum 
behaviour of the integral I{k'^) in eq. ([5]) is dominated by the result of the integration over the 
IR domain, over g2 ^ g,2^ ^nd it can be expanded on 6'^ = p'^/k'^ — 1 around zero with the 
subtraction momentum //2 ^gpt fixed, as explained in ref. [2l|TT], to give 
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*This is the massive gluon propagator where the ghion running mass [31], M(q^), appears to be approximated 
by its frozen value at vanishing momentum, Af (0). It should be also noted that, provided that the gluon propagator 
is to be multiplicatively renormalized, the mass scale, M — M(0), does not depend on renormalization scale, fi^. 



where 
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^ ^ ' (g2 + fe2)2i+2 (g2 + ^2 + ^2)2i+2j (9) 
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Now, the two possible cases for the low-momentum behaviour of the GPDSE solutions will 
be separately analyzed in the following by applying eq. ([8]) to eq. ([1]). 

2.1 Critical case: scaling solution 

When ap < 0, one can perform the change t = q^ /k"^ and integrate over t, the integration in 
eq. ([8]) for any term of eq. ^ being (one by one) convergent as the limit q^/k"^ — >■ cxd is considered 
for the upper bound. Thus, the small-momentum asymptotical behaviour of the r.h.s of eq. (JH) 
is given by 

dt-, TTTTT + (2« + 1) / dt- -r-— 

+ Oid'^), (11) 

while its l.h.s. behaves as: 
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Then, we should conclude both 

(i) that —ap = 1 + ap =^ oip = —1/2, 
(ii) and, given that 



i=0 

one obtains 
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Thus, we know the asymptotic behaviour for the ghost dressing function in this case to be : 

/ in \ 1/2 / 1 \ 1/2 

vr / 10 \ ' / 1 ^ ' 



This is a so-called scaling solution where, in particular, the low-momentum behavior of the 
massive gluon propagator forces the ghost dressing function to diverge at low-momentum through 
the requirement 13 for the power exponent, ap, in (i). 

If ap > is assumed, one would be also left with the contradictory scaling condition that 
ap = — 1 l2l[T7] and should conclude that no solution exists for that case. 



^This is of course a particular case, with aa = 1, of the more general scaling condition: 2olf + olq ~ 0. 



2.2 Regular case: decoupling solution 

When ap = 0, the previous ressumation cannot be done and the integral I{k'^) in eq. ([5]) 
should be consistently expanded in powers of fc^. This case has been deeply studied in ref [TTj . 
where it is found that 
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Then, the ghost dressing function, including its first correction to the leading constant term, 
should behave as |j 
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where 



aT(0) = lim (,2 + M^) ^4^ = m2%!1f1(0)Ah(0), (18) 

g->0 g^ 47r 

such that the eq. (JH) could be satisfied. It should again understood that the subtraction mo- 
mentum for all the renormalization quantities is /i^, j^ eq. (jlSp . a^ = Qr/i'^'^) i^ ^^^ pertur- 
bative strong coupling defined in this Taylor scheme [7], while oIt is the extension of the non- 
perturbative effective charge definition from the gluon propagator [32] to the Taylor ghost-gluon 
coupling [12]. As a consequence of the appropriate amputation of a massive gluon propagator, 
where the gluon mass scale is the same Rl-invariant mass scale appearing in eq. ([6]), this Taylor 
effective charge is frozen at low-momentum and gives a non- vanishing zero-momentum value in 
terms of which the ghost-dressing-function subleading correction can be expressed. 

3 Comparison with numerical results from coupled PT-BFM 
DSE's 



We shall now compare the formulas given by eqs. (|6I17|) with some numerical results for the 
gluon propagator and ghost dressing function. The aim of the comparison is twofold: testing 
the asympotical solution we obtained in the previous section, but also checking the consistency 
of a massive-gluon solution and determining the gluon mass as the best-fit parameter in the 
comparison. In particular, we will consider the solutions of the coupled system of gluon and 
ghost DS equations obtained by applying the pinching technique in the background field method 
(PT-BFM) [13] (see also [33] and references therein) to compare with. This PT-BFM frame- 
work leaves us with an attractive model for gluon and ghost propagators providing quantitative 
description of lattice data |11[33] and giving well account of their main qualitative features: fi- 
nite gluon propagator and finite ghost dressing function at zero-momentum. Futhermore, the 
coupled DSE system can be solved with different boundary conditions (see below), the solutions 



®It should be also noted that eq. H16II imples to take A/^/go "C 1. However, any correction to that approximation 
will not play at the order of the coefficient eqs. (|18|1 . that will keep the same value disregarding that of A/^/go, 
but at the order of the gluon mass, A/^, inside the logarithm, presummably like the UV part of the integral I{k^) 
that should be proportional to k^ and vanish at least like 1/log (go/^QCD)) when go/^QCD S> 1. 



compared with the analytical formula and how their behaviour depends on these boundary 
conditions can be thus properly studied. This last is the main purpose of this note. 

The main feature in the PT-BFM scheme is that the transversality of the gluon self-energy 
is guaranteed order-by-order in the dressed- loop expansion, this leading to a gauge-invariant 
truncation of the gluon DSE [Ijj. In this PT-BFM scheme for the coupled DSE system, the 
ghost propagator DSE is the same as given by eqs. ()ll4p . where the bare ghost-gluon vertex is 
approximated by Hi = 1. The gluon DSE is given by 

X(^2) \9f.u-^\ =q gt,u-qt,qu + i2^{ai)^^ (19) 

where 
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a,= ^...O-.-W. > «4= ,_,^VoV',.,^ • (20) 



In the diagrams of (pO|) for the gluon DSE, eq. p9j) . the external gluons are treated, from the 
point of view of Feynman rules, as background fields (these diagrams should be also properly 
regularized, as explained in (33] )• The last justifies the four field coupling of two background 
gluons and two ghosts leading to the contribution 04. The function 1-l-G defined in ref. [35j can 
be, in virtue of the ghost propagator DSE, connected to the ghost propagator [12]. The coupled 
system is to be solved, by numerical integration, with the two following boundary conditions 
as the only required inputs: the zero-momentum value of the gluon propagator and that of the 
coupling at a given perturbative momentum, /i = 10 GeV, that will be used as the renormaliza- 
tion point. The latter can be done by fixing different values for the boundary conditions, this 
providing us with a family of gluon and ghost propagators solutions so determined. In particu- 
lar, solutions obtained by keeping the zero- momentum value of the gluon propagator fixed (see 
lefthand plots of fig. [1]) while a{p? = 100 GeV ) is ranging from 0.15 to 0.1817 are available |13j 
and can be confronted to the asymptotical expressions derived in the previous section. 

3.1 Decoupling solutions in the PT-BFM scheme 

Then, as the gluon propagator solutions in the PT-BFM scheme result to behave as massive 
ones, the eqs. ()6ll7p must account for the low-momentum behaviour of both gluon propagator 
and ghost dressing function with Hi = 1 and 

aT(0) = aT(//')F|(0)i?(/i2) = ar(/x2)F2(0)M2Ai,(0) , (21) 

«r(A*^) = 5|e(^^)/(47r) being fixed, as a boundary condition, at the moment of the numerical 
integration of the coupled DSE for each particular solution of the family. -B(/i^) and M being 
determined by the best fit of the eq. ([6]) to the numerical solution for the gluon propagator, we 
shall be left with one only free parameter, Fr(0), to account with eq. p7p for the numerical 
solution for the ghost propagator. Furthermore, the zero-momentum values of the ghost dressing 
function, i^/j(0), can be also taken from the numerical integration of the DSE (for any value 
of the a(// = lOGeV)); and these altoghether with the zero-momentum values of the gluon 
propagator, Aj:i(0), and the gluon masses, obtained by the fit of eq. (0) to the numerical DSE 
gluon propatator solutions, provide us with all the ingredients to evaluate, with no unknown 



parameter, eq. (|17p . The gluon masses obtained from the best fits of eq. ([U]) to the numerical 
data (see the left plots in fig. [1]) and that of zero- momentum Taylor effective charge, 0-^(0), 
computed by applying eq. (j2ip . with the zero-momentum ghost dressing function taken from 
numerical data, can be found in tab. [1] 
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Figure 1: Gluon propagators (left) and ghost dressing functions (right) after the numerical integration 
of the coupled DSE system for a(/i = lOGeV) = 0.15, 0.16, 0.17 taken from [T3] . The curves for the best 
fit of eq. © to the gluon propagator data appear as dotted lines in the lefthand plots. In the righthand 
plots, the red dotted lines correspond to apply eq. (|T7)) with the gluon mass obtained from the gluon 
fits and with R = cJt (0)/Af^ determined by the zero- momentum values of gluon propagator and ghost 
dressing function coming from the numerical integration of the DSE system; the same for the black dotted 
lines but retaining only the logarithmic leading term in eq. (J17p by dropping the —11/6 away. 



a(/x) 


ar(0) 


M (GeV) [gluon] 


0.15 


0.24 


0.37 


0.16 


0.30 


0.39 


0.17 


0.41 


0.43 



Table 1: Gluon masses and the zero- momentum non-perturbative effective cliarges, obtained as explained 
in the text, which are applied to describe the gluon propagator and the ghost dressing function numerical 
data with eqs. (|6ll7p . 

Indeed, the expression given by eq. (|17|) can be succesfully applied to describe the solutions 
all over the range of coupling values, a(/i), at /U = 10 GeV. This can be seen, for instance, for 
a = 0.15, 0.16, 0.17, in the right plots of fig. [H where the ghost dressing functions obtained from 
the numerical integration of the DSE's appear plotted with black-solid lines and the evaluation of 
eq. ()17p . as we explained above, with red-dotted ones. Thus, with no parameter to be fitted (Ffi{0) 
is taken from the numerical integration), we nicely reproduce the low-momentum behaviour of 
the ghost dressing function obtained through numerical integration. 

In the right plots of fig. [H black-dotted curves obtained by only retaining the logarithmic 
leading order in eq. (J17p appear also drawn. The big discrepancy they show with respect to the 
numerical integration clearly implies the necessity of the next-to-leading correction in eq. p7|) 
when assessing the gluon mass from the low-momentum ghost propagator. 



3.2 The "critical" limit in the PT-BFM scheme 

There appears to be a critical value of the coupling, Qcrit = oi{^'^) ~ 0.182 with /i = 10 Gev, 
above which the coupled DSE system does not converge any longer to a solution [T^. As a 
matter of the fact, we know from eq. ()14p that the scaling solution implies for the coupling 



2 
9crit 



gW) 



lOvr^ 



3^2(^2)^(^2) ' 



(22) 



where B{^^) is determined by the gluon propagator solution that is supposed to behave as 
eq. ([6]), and A{^^) by the ghost propagator that should behave as 
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where again /i^ is the momentum at the subtraction point. This is shown in ref. pQ, where only 
the ghost propagator DSE is solved there after extracting a gluon propagator from the lattice 
data and applying it to build the kernel of the integral, eq. (I5|), appearing in eq. ([H). In the 
analysis of ref. [1], a ghost dressing function solution diverging at vanishing momentum appears 
to exist and verifies eqs. ()22l23p . while regular or decoupling solutions exist for any a < Ocrit- 

Now, we can perform a more complete analysis by studying again the dressing function 
computed by solving eq. (fT9]) for the different values of the coupling, a = a(^^), at /i^ = 100 
GeV'^ [13]. A ghost dressing function at vanishing momentum, F(0,//^), diverging as q ^^ Ocrit 
had to be expected, one could try the following power behaviour. 
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to describe the vanishing-momentum ghost dressing function in terms of the coupling, a(/.i^). 
The coefficient ^(/i^) should be the positive critical exponent (depending presummably on the 



renormalization point, ^u^) governing the transition from decoupling (a < acrit) to the scaUng 
{a = ttcrit) solutions. 

Our strategy will be to let Ocrit be a free parameter to be fitted by requiring the best linear 
correlation for log[i^(0)] in terms of log[acrit — a]- In doing so, the best correlation coefficient is 
0.9997 for Ocrit = 0.1822, which is pretty close to the critical value of the coupling above which 
the coupled DSE system does not converge any more, and then we obtain 
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Figure 2: (a) Log- log plot of the zero- momentum values of the ghost dressing function, obtained by 
the numcrieal integration of the coupled DSE system in the PT-BFM scheme, in terms of ctcrit ~ ct- 
a = a{^ = fOGeV), the value of the coupling at the renormalization momentum, is an initial condition 
for the integration; while Ocrit is fixed to be 0.1822, as explained in the text, by requiring the best linear 
correlation, (b) Gluon propagator solutions in terms of q^ for the same coupled DSE system for different 
values of a{iJ. = lOGeV), all very close to the critical value, ranging from 0.18 to 0.1817 . 

In fig. [5J(a), the log-log plot of Fr{0) in terms of Ocrit — a is shown for Ocrit = 0.1822, where 
the linear behaviour corresponding to the best correlation coefficient can be strikingly seen 
and the negative slope indicates a zero-momentum ghost propagator diverging as a — >■ Ocrit- 
However, no critical or scaling solution of the coupled DSE system seems to appear with a 
massive gluon propagator as solution of the coupled DSE system in the PT-BFM, although the 
decoupling solutions obtained for any a < ctcrit = 0.1822 seem to approach the behaviour of a 
scaling one when a — > acrit- The absence of the scaling solution can be well understood by 
analysing eq. ()19p . As explained in ^, after the appropriate regularization and renormalization, 
the contribution of 03 to the inverse of the gluon propagator, its momentum vanishing, will be 
dominated byL| 
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'\{q^) {F\q^) - I) 



(26) 



where qq is again some UV cut-off above which the ghost dressing function can be taken to be 
perturbative. Provided that one deals with a decoupling solution, the ghost dressing function 



^The regularisation procedure in [J implies the subtraction of the perturbative part, as well as we evaluate 
eq. (m for two different scales and subtract them in order not to have to deal with any UV cut-off in eq. ^. 
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reaching some constant as g^ — )> 0, this contribution is finite and neghgible (the same happens for 
04). On the other hand, had we considered the scaling solution, the ghost dressing function would 
behave as Xjq and would lead to a divergent contribution and to a vanishing gluon propagator 
at vanishing momentum. Then, a massive gluon propagator, as lattice solutions points to and 
as required in ref. [^113]. cannot appear as a scaling solution. The same has been already proven 
in literature by applying different truncation schemes and also on general grounds, and indeed 
agrees with the very well-known infrared behaviour obtained from the coupled DSE system in 
refs. [316], where an unique scaling solution with ap ~ —0.595 and ag = —2oip ~ 1.190 emerges, 
the gluon propagator vanishing thus at zero-momentum, as the inverse of the zero-momentum 
ghost dressing is assumed to vanish. However, for the sake of completeness, a very general 
argumentation addressing this issue is presented in appendix [XI We should recall at this point 
that the numerical analysis of the ghost propagator DSE in ref. [T] left us with a divergent 
ghost dressing function for the critical value of the coupling, even after assuming a finite gluon 
propagator at zero momentum. However, this resulted from a partial analysis where we dealt 
only with the ghost propagator DSE and not with the gluon propagator one. Thus, we did not 
solve the gluon DSE and taking a massive gluon propagator from the lattice to build the ghost- 
DSE kernel does not prevent from obtaining a ^^wrong" scaling solution indeed not satisfying the 
coupled DSE system. 

When approaching the critical value of the coupling, the gluon propagators obtained from 
the coupled DSE system in PT-BFM must be also thought to obey the same critical behaviour 
pattern as the ghost propagator. In the PT-BFM, the value at zero-momentum being fixed by 
construction [^[[TS]. one should expect that, instead of decreasing, the gluon propagator obtained 
for couplings near to the critical value increases for low momenta: the more one approaches the 
critical coupling the more it has to increase. This is indeed the case, as can be seen in fig. [21^b). 
This implies that, near the critical value, the low momentum propagator does not obey eq. ([6|) 
and that consequently eq. ()17p does not work any longer to describe the low momentum ghost 
propagatoilj. 

Finally, one can pay attention to the critical value of the coupling, Ocrit = 0.1822, and try to 
make a comparison with the physical strong coupling values in order get some idea of whether the 
current data can exclude or not this critical behaviour. Although the experimental PDG world 
average of the strong coupling in the MS scheme, a-^{Mz) = 0.1184(7) [36], can be propagated 
from the Z^ boson mass down to ^ = 10 GeV to give aj^(10 GeV) = 0.179(2), that incidentally 
lies on the right ballpark of the above critical value, such a comparison is meaningless because 
our coupling corresponds to one in MOM Taylor-scheme for zero number of flavours. One can 
use instead the available perturbative four-loop formula describing the running of the coupling 
in Taylor-scheme to estimate Aqcd in this particular scheme, then perform the conversion to 
MS (see for instance eqs. (22,23) of the first reference in [7]) and thus obtain the value quoted in 
tab. 13. 2i Of course, it would be again meaningless to compare this last value with the one for 
Aj^ that can be obtained from the PDG value for aj;jg(M^), also quoted in tab. 13.21 but we 
can refer the comparison to the lattice Yang-Mills determinations of the same parametecj, as for 
instance the two of them included in tab. 13.21 Thus, the lattice estimates of Aj^ appear to lie 
clearly below this critical limit for the PT-BFM DSE in pure Yang-Mills. However, as no quark 
flavour loops effect have been incorporated in our DSE, eq. (|19p . we cannot neither compare 
with the physical strong coupling nor conclude whether the critical limit can be allowed in the 



*This is only true for the next-to- leading contribution of eq. {TTJ, the leading one being only determined by 
the zero-momentum gluon propagator still works. 

''it should be noted that the procedures for the lattice determination of A^,^ mainly work in the UV domain, 
where IR sources of uncertainties as the Gribov ambiguity or volume effects are indeed negligible. In fact, there 
are unquenched lattice determinations with Nf = 5 staggered fermions for the strong coupling [33 which are 
pretty consistent with the PDG value. 
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^^real world''. 












MS.crit 


Air" m 


^T m 


A^-'m 




434 MeV 


238(19) MeV 


244(8) MeV 


213(9) MeV 



Table 2: The critical value of Aj^ in pure Yang-Mills inferred from acrit — 0.1822 (first column), lattice 
estimates for Yang-Mills A^^jg taken from literature (second and third eolumns) and the one obtained 
from the PDG value of a^^{Mz) by applying a four- loop perturbative formula for the running of aj^g 
with iV/ = 5. 



4 Conclusions 



The ghost propagator DSE, with the only assumption of taking Hi{q, k) from the ghost-gluon 
vertex in eq. ([3]) to be constant in the infrared domain of q, can be exploited to look into the low- 
momentum behaviour of the ghost propagator. The two classes of solutions named "decoupling" 
and "scaling" can be indentified and shown to depend on whether the ghost dressing function 
achieves a finite non-zero constant {ap = 0) at vanishing momentum or not {ap 7^ 0). The 
solutions appear to be dialed by the size of the coupling at the renormalization momentum 
which plays the role of a boundary condition for the DSE integration. The low-momentum 
behaviour of the decoupling solutions results to be regulated by the gluon propagator mass and 
by a regularization-independent dimensionless quantity that appears to be the effective charge 
defined from the Taylor-scheme ghost-gluon vertex at zero momentum. 

In this note, we have studied the solutions of coupled ghost and gluon propagator DSE 
in the PT-BFM scheme and demonstrated that the asymptotic decoupling formula {ap = 0) 
successfully describes the low-momentum ghost propagator. The model applied for the massive 
gluon propagator is also verified to give properly account of the gluon solution, at least for 
momenta below 1 GeV (and for a coupling not very close to the critical point). Although we 
argued that a massive gluon propagator implies that the ghost dressing function takes a non- 
zero finite value at vanishing momentum, we also show that the zero-momentum ghost dressing 
function tends to diverge when the value of the coupling dialing the solutions approaches some 
critical value. Such a divergent behaviour at the critical coupling seems to be the expected one 
for a scaling solution (where, if the gluon is massive, ap = —1/2). If we consider the zero- 
momentum value of the ghost dressing function as some sort of ^^ order parameter'^ indicating 
whether the ghost propagator low-momentum behaviour is suppressed {ap = and finite ghost 
dressing function) or it is enhanced {ap < and divergent ghost dressing function) , the strength 
of the coupling computed at some renormalization point seems to control some sort of transition 
from the suppressed to the enhanced phases for the ghost propagator DSE solutions in the 
PT-BFM scheme. The last only takes place as some critical value of the coupling is reached. 
Neverteless, it can be proven that, as far as the gluon is massive, the scaling behaviour for the 
Yang-Mills propagators appear not to be a solution but an unattainable limiting case for the 
PT-BFM DSE solutions. 

Acknowledgements: The author is particularly indebted to Ph. Boucaud, J.P Leroy, 
A. Le Yaouanc, J. Micheli and O. Pene for very fruitful discussions at the initial stages of 
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the work and to J. Papavassiliou and A.C. Aguilar also for very valuable discussions and com- 
ments, and specially for providing me with some unpublished results which were exploited in 
this paper. J. R-Q also acknowledges the Spanish MICINN for the support by the research 
project FPA2009-10773 and "Junta de Andalucia" by P07FQM02962. 

A No scaling solution with massive gluons 

We consider the conventional gluon self-energy, U.^'^{q), contributing to the gluon DSE: 



.^'^^. 



^f^u{q) - i: ssAiSiQOQSi^SiSh +9 . o o qV^s-^q o . + ^:ia0- - -O- - - W^ 



2 viiiA.jioiiii_iioiii,»».iisi^ 2 oi_2_fi_ai:S3i_fi_SLJL^ 




,-o~. 




(27) 



where the yellow bullets stand for full vertices and propagators. 

After only assuming that the ghost-gluon vertex form factor Hi is constant (the Taylor non- 
renormalization theorem tells us that the bare ghost-gluon vertex is finite) , We showed in section 
[2] that a massive gluon propagator unequivocally implies a ghost dressing function diverging as 
1/q at vanishing momentum. Then, the ghost-loop contribution to the gluon self-energy, eq. ([2^ 
with vanishing external momentum, k, is dominated by 



9, 



T n.^ ., X ^ •> , / d'l „2 f. (k ■ qr\ F{q') F{{q - kf) 



,Ak)x ^...o-W - l^^'^'i'-^-^) .2 



g2 {q - kf 



3/2 



go „2 



' (g2 + A;2)3/2 

where go is the momentum scale which the ghost dressing function is assumed to take its asymp- 
totic infrared form below. For not to have to deal again with any UV regularization cut-off, we 
can consider the two momenta p, k such that k^ <^ p'^ < q^ and subtract the gluon DSE for these 
two momenta. Then, when A;^ ^ while p'^ is kept fixed, one would have for the transversal 
gluon propagator 

go „2 

dq , , \, .,, , (29) 



AR(fc2) Anip^) 7o "(g2 + fe2)3/2 

where we only account for the dominant part of the ghost-loop contribution. This contribution 
diverges thus logarithmically with a constant ghost-gluon vertex again as the only assumption 
made, while computing the contributions coming from other diagrams would leave us with the 
neccessity to make some new assumption about the full gluon vertices. 

Thus, to avoid a diverging gluon self-energy, we need to invoke new contributions from the 
other diagrams in eq. ()27p . also diverging logarithmically as the external momentum vanishes, to 
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cancel that from eq. (|28p . Otherwise, such a divergent behaviour of the gluon self-energy would 
lead the inverse of the gluon propagator to diverge, and the gluon propagator consequently to 
vanish, at zero-momentum. 
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